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Abstract 

We  study  the  nature  of  triplet  states  in  correlated,  quasi  one  dimensional  bands, 
with  emphasis  on  the  half-filled  case  relevant  to  conducting  polymers.  To  incorpo¬ 
rate  both  electron-phonon  and  electron-electron  interaction  effects,  we  use  a  Peierls- 
extended  Hubbard  Hamiltonian,  which  we  solve  for  finite  systems  via  (numerically) 
exact  diagonalization  (Lanczos)  techniques.  We  extend  our  results  to  polymers  using 
both  standard  finite-size  extrapolation  techniques  and  a  novel  boundary  condition 
averaging  scheme.  First,  we  examine  the  nature  of  the  lowest  triplet  state,  focusing 
on  the  crossover  from  the  free  electron  picture  of  neutral  soliton  pairs  to  the  strongly 
correlated  limit  of  spin  excitations.  We  then  examine  the  optical  absorption  spectra 
corresponding  to  (allowed)  triplet-triplet  transitions;  for  this  purpose  it  is  important 
that  our  novel  boundary  condition  averaging  method  maintains  sharp  gap-edge  fea¬ 
tures  while  reducing  spurious  rapid  frequency  variations  p'-oduced  by  the  finite-size 
approximations  to  the  band-to-band  continuum.  We  compare  our  results  to  data 
on  triplet  transitions  in  finite  polyenes  and  discuss  possible  additional  experimental 
implications,  including  the  interpretation  of  recent  spin-dependent  photomodula¬ 
tion  (SDPM)  experiments  in  both  trans-{CH)x  and  non-degenerate  ground  state 
conjugated  polymers. 


I.  Introduction 


No  interpretation  of  the  excited  state  spectra  and  in  particular  of  the  dynamics  of  photoex¬ 
citations  in  finite  polyenes  and  conducting  polymers  can  be  considered  complete  unless 
it  correctly  captures  the  important  role  of  the  triplet  stales.  In  (derivatives  of)  finite 
polyenes,  triplet  states  are  known  to  be  of  considerable  importance  in  the  processes  of 
vision  [1]  and  of  photosynthesis  [2],  in  that  they  control  crucial  relaxation  paths  after  elec¬ 
tronic  excitation.  In  conducting  polymers,  triplets  appear  essential  to  understanding  the 
nature  of  long-lived  excited  states  and  specifically  of  the  spectra  observed  in  photoinduced 
photoabsorption.  From  a  theoretical  perspective,  knowledge  of  triplet  states  and  their 
properties  can  provide  additional  insight  into  the  relative  importance  of  electron-electron 
(e-e)  and  electron-phonon  (e-p)  interactions  and  can  help  determine  the  parameters  in  the 
theoretical  models  of  these  materials. 

Fortunately,  a  number  of  recent  experimental  results  have  begun  to  clarify  the  prop¬ 
erties  of  triplet  states,  including  the  energy  levels  and  their  (linear)  optical  absorption 
properties.  In  finite  polyenes,  earlier  work  using  pulse  radiolysis  and  flash  photolysis 
[1,  2]  to  observe  the  very  weak  (forbidden)  singlet-triplet  in  addition  to  the  (allowed) 
triplet-triplet  transitions  has  be  supplemented  by  the  results  of  ultra  high  resolution  spec¬ 
troscopy  [3]  and  of  embedding  the  polyenes  in  “  spectroscopic  matrices”  [4]  which  enhance 
the  forbidden  transitions.  In  conducting  polymers,  novel  experimental  techniques  includ¬ 
ing  “optically  detected  magnetic  resonance”  (ODMR)  [5,  6,  7,  8],  photcinduced  ESR  [9], 
and  “spin-dependent  photomodulation”  (SDPM)  [10,  11]  have  begun  to  produce  definitive 
information  on  the  location'  and  properties  of  triplet  excitations. 

Theoretical  studies  of  triplets  in  these  systems  have  to  date  been  less  extensive. 
Within  the  framework  of  the  SSH  model  [12],  Su  [13]  showed  that  the  lowest  triplet 
excitation  was  expected  to  be  a  pair  of  neutral  solitons  with  aligned  spins.  Hartree- 
Fock  [14],  renormalization  group  (15,  16]  and  configuration  interaction  [17]  studies  have 
suggested  that  in  models  in  which  electron-electron  interactions  are  present,  triplets  could 
play  an  especially  important  role  as  “constituents”  in  the  covalent  singlet  excitations,  such 
as  the  “cplebrat(xJ”  2' Ag  state. 

However,  these  intriguing  pioneering  studies  have  not  been  follow'ed  by  definitive 
analyses  which  .simultaneously  (1)  incorporate  completely  the  full  effects  of  the  possible 
geometric  relaxations  of  triplets  and  (2)  solve  exactly  the  true  many-body  problem  that 
arisers  wIkui  <>  e  interactions  are  properly  included.  Our  goal  in  the  present  article  is  to 
licgin  such  a  definitive  analysis  by  studying  via  an  exact  many-body  method  a  model 
which  iiicorpoiates  both  e-p  and  e-e  interactions  of  arbitrary  strcuigths.  We  begin  in 
Section  11  by  introducing  the  Peierls-extended  Hubbard  mode]  m  the  form  appropriate 
/irm.s  p.ilvai  ftylene  { trans-{(l  ll)x:)  ‘"md  discus.sing  the  numerical  exact  diagunalization 
(Lanezbs)  tei  hnicpie  wc-  use  for  solving  the  model.  We  also  explain  our  “boundary  condi 
tion  averaging  method  for  e.xtracting  estimates  ol  energies  ai..l  otl.cr  propcties,  wliich 
estimati';.  an-  i  'inside! ai)l\-  less  sensitive  to  finite  size  effi'cts  tlian  (dios<’  obtained  by  i  on 
ventioiial  iii'-arr.  In  Sintiou  Ill  we  review  tin*  analytiiallv  tra.c  table  weak  {SSH  ninile!) 
and  ilioii!'  i  li'Mse-  i  ■  ;  -  ::,,,.lc;;  •o.qil.ug  limii-'  <md  t.i.si  u.-^s  I  in'  exjiectat  ions  tor  tiiplel 
piopoi  I  K';.  li;  '  in'se  limil  mg  va-'es.  Se(  tion  |\'  contains  tiie  dolaiF  ol  om  iinmcrica!  ii  :  alts 


for  the  energies,  geometries,  and  optical  absorptions  of  both  triplet  and  (for  comparative 
purposes)  singlet  states.  In  Section  V  we  examine  the  possible  role  of  triplet  states  in 
recent  spin-dependent  photomodulation  experiments  and  the  relation  of  the  triplet  states 
to  the  (optically  forbidden)  2^  Ag  state.  Finally,  in  Section  VI  we  discuss  a  number  of 
unresolved  issues  and  future  challenges  for  both  theory  and  experiment. 


II.  The  Peierls'Hubbard  Model  and  Lanczos  Method 


In  the  context  of  conducting  polymers,  the  one-dimensional  Peierls-extended  Hubbard 
Hamiltonian  [18]  provides  a  theoretical  framework  capable  of  treating  both  electron- 
phonon  (e-p)  and  electron-electron  (e-e)  interactions  of  arbitrary  strengths.  In  the  specific 
case  of  frans-polyacetylene  {irans-{C H)^),  the  Hamiltonian  is 

e  t  ^  t  t  t 

where  c\^{cea)  creates  (annihilates)  an  electron  in  the  Wannier  orbital  at  site  f,  nt  =  -|- 

nei,  where  nt„  =  is  the  density  operator,  and  „ct+\  a  +  c\^x  gCt  „)  is 

the  electron  hopping  operator,  though  it  may  also  usefully  be  view'ed  as  a  "bond-charge” 
operator  [19|.  St  represents  the  projection  along  the  chain  axis  of  the  relative  displacement 
between  the  {CH)  units  at  sites  £  and  £-f- 1  and  pt  the  momentum  of  the  (CH)  unit  at  site 
i  and  mass  M.  The  single  particle  parameters  of  the  model  are  the  hopping  integral  to  for 
the  uniform  (CH)  ionic  lattice,  the  electron-phonon  coupling  o  describing  the  modification 
of  the  hopping  between  adjacent  sites  due  to  the  distortion  of  the  underlying  discrete 
lattice,  and  the  lattice  spring  constant  K  representing  the  cost  of  distorting  the  lattice  of 
[C H)x  moieties  due  to  interactions  not  otherwise  explicitly  included  in  the  model,  such 
as  the  core  electron  repulsion.  Since  in  some  of  our  calculations  we  shall  be  interested 
in  the  comparative  geometries  of  the  ground  and  excited  states,  rather  that  using  the 
usual  constraint  =  0,  we  have  defined  oq  to  be  the  value  needed  in  order  that  in 

the  uniformly  dimerized  ground  state  the  sum  is  zero.  Using  this  value  of  oq  allows  to  to 
determine  the  excited  geometries  in  a  manner  consistent  with  that  in  the  ground  state. 
As  in  almost  all  studies  of  this  model,  we  work  in  the  adiabatic  limit,  ignoring  the  lattice 
momentum  term  in  Eq.  (1)  and  treating  the  lattice  coordinates  as  classical  variables  whose 
values  are  determined  by  minimizing  the  electronic  energy  self-consistently. 

The  electronic  many-body  terms  due  to  Coulomb  repulsion  are  parameterized  by  the 
conventional  Hubbard  IJ  and  V^’,  de.scribing  the  on-site  and  nearest-neighbor  interactions, 
r(*spectively.  In  the  limit  IJ  --  V''  =  0  (dominant  e-p  interactions).  Piq.  (1)  reduces  to  the 
familiar  Su-SchriefTer-Heeger  (SSll)  model  [12]  of  lrans-[C U)x-  In  the  a  =  0  limit  (no  e-p 
but  dominant  e-e  interactions),  it  reduces  to  the  extended  Hubbard  model  [19,  20). 

Considerable  experimental  evidence  (reviewed  in  [18])  suggests  that  parameter  val- 
ue.s  appropriate  tor  (  ondiirtiiig  polymers  lie  in  the  intermediate  coupling  regime:  C'  Ain. 
wIhmi-  Ain  is  the  tu!!  omi;h:  jau  in  ie  bandwidUi.  To  obtain  definitive  results  in  this  pa 
rameter  regime,  one  must  use  evac  I  many  body  techuiqui's.  We  ado[>l  adopt  the  hanrzbs 
c.x.u  t  diagonaiization  to  hmcjii''.  which  we  describe  very  briefly  below  bei  aiisc"  of  some  sig 


nificant  modifications  which  simplify  considerably  the  meaisurement  of  optical  absorption 
spectra. 


For  studies  of  the  optical  transitions,  we  shall  need  the  current  operator 

jt.f+1  =  i(<o  -  aS()  a)  .  (2) 

<T 

The  Fourier  transform  of  is 

=  .  (3) 

The  optical-absorption  coefficient  a(ui)  is  given  by 

a(")  =  -  E  l(”>|J',|0)|’<(u,  -  (£„  -  £„))  .  (4) 

^  m 


In  our  variant  of  the  Lanczos  algorithm,  we  start  generating  the  basis  set  by  normaliz¬ 
ing  some  trial  wavefunction,  which  we  typically  express  in  terms  of  real  space  occupations. 
We  have  usually  chosen  to  start  from  a  random  wavefunction,  subject  only  to  the  con¬ 
straint  that  5^  =  0,  where  St  is  the  z-component  of  the  total  spin.  Since  the  Hamiltonian 
contains  no  (net)  spin  flip  operators,  this  constraint  is  preserved  throughout  the  calcu¬ 
lations.  Choosing  a  random  starting  wavefunction  does  not  cost  much  extra  computing 
time  in  most  cases  and,  importantly,  statistically  prevents  one  from  picking  a  wavefunction 
with  the  wrong  symmetry. 

The  Lanczos  procedure  is  to  generate  the  matrix  elements  of  the  Hamiltonian  in  a 
basis  that  is  built  up  from  this  trial  state.  The  basis  is  incremented  one  wavevector  at  a 
time  oy  operating  on  the  leist  basis  state  with  the  Hamiltonian  and  then  orthonormalizing 
the  product  to  all  previous  basis  states.  Notice  that  by  construction  {j\H]i)  is  tridiagonal 
[21,  22].  Within  this  basis,  that  linear  combination  of  states  with  the  lowest  expected 
energy  forms  the  estimate  of  the  ground  state.  The  basis  will  probably  not  be  complete: 
we  generally  will  truncate  it  when  we  have  run  out  of  computer  memory,  when  the  estimate 
of  the  ground-state  energy  stops  dropping,  or  when  the  residual  of  a  new  state,  after 
orthogonalizing  to  other  components,  is  negligible. 

Importantly,  the  calculations  of  optical  absorption  spectra  - —  and  hence  also  gaps, 
conductivities,  and  susceptibilities  —  are  easily  carried  out  by  generating  H  in  a  suitable 
basis.  For  the  optical  absorption,  we  first  calculate  the  ground  state,  |v!^o)  using  the  tech¬ 
nique  described  above.  We  then  generate  a  new  basis,  using  Jg\d>o)  as  the  first  state  and 
obtaining  subsequent  states  by  the  application  of  //.  Once  the  b2isis  has  been  generated, 
the  problem  of  calculating  a(u;)  reduces  to  finding  the  spectral  weight  of  this  first  state 
for  the  tridiagonal  Hamiltonian.  In  essence,  this  amounts  to  determining  the  spectrum 
by  measuring  the  moments  of  //'*  and  using  the  curnulant  expansion. 

The  final  theoretical  ingredient  in  our  analysis  is  somewhat  technical  and  involves 
the  use  of  a  novel  “boundary  condition  averaging”  technique  for  effectively  reducing  thi’ 
finite  size  corrections  to  optical  gaps  and  spectra  anri  thereby  enhancing  tlu'  conviugence 


•t 


of  finite  system  results  to  those  of  the  infinite  system.  Readers  interested  in  the  technical 
details  should  consult  the  original  references  (22,  23].  Here  we  note  just  that  by  “boundary 
condition  averaging”  we  mean,  for  example,  that  the  total  energy  of  the  system  is  viewed 
as  the  average  of  the  energies  derived  separately  for  each  of  several  boundary  conditions; 

^  =  7!-  E  I  (5) 

where  N^fc  is  the  number  of  boundary  conditions  used  and  Xf,c  is  a  normalized  weighting 
factor,  which  may  depend  on  which  quantity  we  are  studying.  Similarly,  the  optical 
absorption  is  viewed  as  the  average  of  absorptions  from  the  separate  boundary  conditions. 

The  value  of  the  self-consistent  adiabatic  lattice  distortion  is  found  by  minimizing 
the  total  energy 

dEldSt  =  0  .  (6) 

At  fixed  /f,  to  find  the  grou-;d  state  lattice  distortion  we  use  the  Lanczos  procedure  iter¬ 
atively  to  calculate  the  new  guess  at  the  minimum  energy  distortion  from  the  expectation 
value  of  the  bond  chaxge:  this  is  just  the  familiar  “self-consistency  condition” 

K{St  —  B.)  —  —a{Bf(^.i)  (7) 

with  the  proviso  that  it  be  appropriately  averaged  over  the  boundary  conditions. 


III.  Analytic  Limits 


We  can  gain  considerable  insight  into  our  later  numerical  results  -  as  well  as  relate  to 
earlier  work  in  the  literature  -  by  considering  the  two  essentially  analytic  limits  of  weak 
and  strong  coupling. 


A.  The  Weak  Coupling  Limit 

In  accordance  with  the  common  usage,  by  the  “weak  coupling”  limit  of  the  Peierls- 
extended  Hubbard  model  we  mean  the  caise  in  which  (e-e)  interactions  are  weak,  so  that 
an  independent  particle  mode!  provides  an  accurate  approximation.  For  U  —  V  —  0, 
the  Hamiltonian  in  Eq.  (1)  reduces  to  the  SSH  model  [12],  for  which  the  solitonic  nature 
of  the  triplet  slate  has  been  established  by  Su  (I3j.  We  summarize  briefly  the  results. 
Working  in  momentum  space,  one  obtains  the  grouiui  states  of  the  singlet  and  triplet  by 
filling  the  single  particle  (band  theory)  energy  levels  up  to  the  Fermi  surface.  The  singlet 
ground  state  corresponds  to  the  occupancy  show-ii  in  Fig.  la,  m  which  all  levels  below 
the  (Peieris)  band  gap  are  filled.  Momentarily  ignoring  relaxation  of  the  lattice  in  the 
triplet  slate,  we  sec*  that  the  lowa*st  triplet  would  be  as  sketched  in  F'lg.  lb,  since  two 
electrons  of  th<’  same  spin  cannot  share  a  single  band  stale.  Hence  m  this  ap[)roxiination 
the  energy  of  th<‘  lowest  trijilet  would  be  /v”  -s  2A|),  so  that  t.he  triplet  would  occur  <il 
th('  same  energy  a.s  I  he  o[)l  ic.il  gap  llowevac,  a.s  shown  by  Sii  { 1 .1],  vvmI  bin  the  S.SH  model, 
the  trm'  lowest  l,rij)let  stale  in\o|ves  a  rela.xation  of  the  lattice  to  !orm  a  jiair  of  (widely 


separated)  neutral  solitons  (which  is  allowed  topologically).  Hence  the  schematic  single 
particle  energy  levels  are  as  shown  in  Fig.  2.  The  energy  of  this  relaxed  triplet,  which  is 
the  true  triplet  ground  state  in  the  weak-coupling  limit,  is 


Ej-  —  E$q^  — 


(4Ao) 


Thus  the  triplet  lies  below  the  optically  allowed  (*jB„)  state.  Further,  the  geometry  of 
the  relaxed  triplet  state  involves  a  soliton-cintisoliton  lattice  distortion.  Finally,  the  lowest 
triplet-triplet  optical  transition  will  involve  promotion  of  one  of  the  occupied  midgap 
states  into  the  conduction  band  eind  will  thus  occur  at  At  =  Aq  =  1/2 A5,  where  we  have 
denoted  the  singlet-singlet  optical  gap  by  A5. 

For  later  purposes,  we  note  two  specific  features  of  this  weak  coupling  limit.  First, 
the  role  of  “Jahn-Teller”  and  “anti-Jahn-Teller”  systems  is  reversed  between  the  singlet 
and  triplet  states:  that  is,  the  triplet  states  (in  the  absence  of  a  gap)  are  degenerate  at  the 
Fermi  surface  for  4N  -f  2  systems  with  periodic  (or  4  N  anti-periodic)  boundary  conditions. 
This  can  affect  the  interpretation  of  numerical  results  on  small-sized  systems,  at  least  for 
weak  coupling.  Second,  for  (forbidden)  optical  transitions  between  the  singlet  ground  state 
and  the  triplet  state  -  such  as  would  be  visible  in  the  high  resolution  spectroscopy  [3]  or  in 
the  “spectroscopic  matrices”  [4]  -  one  would  expect  to  observe  the  unrelaxed  form  of  the 
triplet  state  because  of  the  essentially  instantaneous  nature  of  the  transition.  However, 
for  allowed  triplet-triplet  transitions  -  as  seen  in  any  of  the  photoinduced  measurements 
-  the  relaxed  triplet  state  should  be  play  the  role  of  the  initial  state. 
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l-'igurc  I  Caption:  Schematic  representations  of  the  momentum  .spat  r  orbital  o<  i  n 
paiiry  for  tlie  unrelaxcd  (perfectly  dimerized)  (a)  singlet  and  (b)  triplet  gronml  sl.at.-s 
'I’his  .situation  i.s  appropriate  for  the  “weak  coupling"  {U/tn  -  0)  limit 


n 

n 


Figure  2  caption:  Schematic  representation  of  the  momentum  space  orbitaj  occu¬ 
pancy  of  the  triplet  ground  state  after  relaxation  of  the  lattice;  the  “two-soliton”  lattice 
deformation  leads  to  two  mid-gap  states,  which  are  occupied  as  for  two  neutral  solitons. 
This  situation  is  appropriate  for  the  “weak  coupling”  ((//to  =  0)  limit. 


B.  The  Strong  Coupling  Limit 

We  consider  next  the  “strong  coupling  limit”,  by  which  we  mean  that  the  Coulomb  inter¬ 
action  terms  are  much  larger  than  all  single  particle  energies.  Working  for  simplicity  with 
C  =  0,  this  implies  that  (//to  — +  oo.  Since  in  this  limit  there  can  be  no  doubly-occupied 
“sites”  [no  =  0),  it  is  most  natural  to  work  in  “real”  (ic,  coordinate)  space  and  to  use 
valence  bond  diagrams  to  describe  the  states.  As  is  well  known,  in  the  absence  of  dimer 
ization  (a  =  0)  the  Hubbard  model  at  half-filling  in  this  limit  becomes  equivalent  to  an 
antijerromagnetic  Heisenberg  spin  system,  with  effective  coupling  Jrjj  ~  IqIH  Again  we 
.summarize  briefly,  using  two  figures,  some  familiar  results  (sec,  eg,  []8j  for  a  more  di’tailed 
survey  o!  results  in  this  limit.) 


In  Fig.  3  we  show  schematic  real  space  electron  distributions  for  the  singlet  and 
triplet  ground  states  in  the  strong  coupling  limit.  In  the  singlet  slate,  loosely  six-aking 
all  the  spins  are  aligned  antiferromagnetically,  whereas  in  the  triplet  slate  there  are  two 
bonds  which  are  not  antiferromagnetic:  hence  the  energy  of  the  triplet  state  in  this  limit 
is  estimated  to  be  Ey  ~  ~  tl/U  0.  Strictly  speaking,  the  Neel  state  .shown  in 

Fig.  3  a  “cartoon”,  since  the  actual  ground  state  can  have  only  local  (not  long-range)  spin 
order.  However,  for  purposes  of  estimating  energies,  this  naive  picture  is  adequate. 
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F'igure  3  caption:  Schematic  representation  of  the  distribution  of  electrons  in  real 
space  in  the  “strong  coupling”  [U jlo  oo)  limit.  The  Neel  state  shown  is  a  cartoon;  the 
actual  ground  state  hais  only  local  (not  long-range)  order. 

Figure  3  suggests  that  in  the  strong  coupling  limit  the  triplet  state  is  a  single  localized 
excitation,  in  contrast  to  the  weak  coupling  triplet,  which  is  composed  of  two  solitons. 
However,  as  shown  schematically  in  Fig.  4,  the  particular  real  space  configuration  with 
three  parallel  spin.s  adjacent  to  <’acli  other  is  in  fact  degenerate  with  (an  infinite  inimber 
of)  other  configurations  in  which  two  pairs  of  two  parallel  spins  are  separated  by  an 
arbitrary  di.stance.  Since  thc.se  configurations  arc  degenerate  to  leading  order,  the  actual 
ground  state  will  be  a  mixture  of  them.  In  the  terminology  of  the  Heisenberg  spin  system, 
the  triplet  .state  consists  actually  of  tv)o  independent  “domain  walls”,  which  mu.sl  each 
lhercfoi(>  have  spin  1/2.  Note  also  that  the  “sense”  of  the  antiferromagnetic  ordering 
i<‘,  th<‘  :i:(--l)"  changes  as  one  pa.s.ses  through  one  of  the  domain  walls.  In  short,  the 
.'-trong  ( oiijding  trijrlet  (in  the  limil  of  zero  dimerization)  also  consists  of  two  iiKiepenrlent 
e.o  itations  anri  wlnm  the  e  p  coupling  is  restored,  we  expert  the  geometry  to  refiiM  t  this. 
.A:-  iK/’cd  .ibove,  tile  energy  of  the  lowest  triplet  approaches  zero,  /v/-  ~  injl  ‘  0 


S 


Further,  the  the  triplet-triplet  optical  gap  can  readily  be  deduced  by  noting  that  the 
matrix  element  of  the  current  operator  (recall  Eqns.  (2)-(4))  in  for  both  singlet-singlet 
and  triplet-triplet  transitions  connects  subspaces  with  rtjp  —  0  to  those  with  no  =  1  and 
thus  involves  an  energy  gap  ~  [/;  hence  Ay  =  As  ~  [/  in  strong  coupling. 


(e» 


Figure  4  caption;  Schematic  representation  of  some  of  the  degenerate  real  space 
electron  distributions  contributing  to  the  triplet  state  in  the  “strong  coupling”  ((//to  — ♦ 
oo)  limit.  Note  that  each  has  (tc  leading  order)  energy  Et  =  m 

To  explore  the  intermediate  coupling  regime  and  the  crossover  from  weak  to  strong 
coupling,  we  turn  in  the  following  section  to  our  Lanczos  approach. 

IV.  Numerical  Re.sults 

I'o  discuss  our  numerical  results  we  begin  by  defining  in  Fig.  5  the  various  observables 
we  will  study.  As  shown  in  that  figure,  these  include  singlet  and  triplet  energies,  opti¬ 
cal  gaps,  and  optical  absorption  spectra  as  functions  of  photon  frequency.  To  organiz'- 
our  presentation,  we  divide  it  into  three  subsections  covering  energetics,  geometries,  and 
optical  absorption. 
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Figure  5  caption:  Schematic  diagram  of  the  various  energy  levels  and  transition 
energies  discussed  in  the  text.  Energy  levels  shown  include  the  ground  state  energy 
{Es  ~  E{VAg)),  the  first  allowed  singlet  optical  state  energy  (£(’J9u)),  the  (relaxed 
lattice)  triplet  ground  state  energy  {E^  =  E{l^B^)),  and  the  “second  A,”  state  en¬ 
ergy  (£(2Mg)).  Transitions  shown  are  the  (allowed)  singlet-singlet  gap  (As)  and  optical 
absorption  coefficient  (as(a;))  and  the  (allowed)  triplet-triplet  gap  (Ar)  and  optical  ab¬ 
sorption  coefficient  (a7-(tci)). 


A.  Energetics 

In  Fig.  6  we  plot  the  electronic  contributions  to  the  energies  of  the  ground  stale  and 
several  important  excited  stale.s  as  functions  of  the  Hubbard  IJ  parameter;  the  states 
shown  include  the  lowest  triplet  (l^i?„),  the  second  singlet  “A,”  .state  (2‘/lg),  and  the 
“optical  state"  (2*  B^).  The  figure  is  plotted  for  an  eight  site  system  with  open  iioundary 
conditions,  corresponding  to  (idealized)  octatetraene.  Note  that  the  minimum  value  of 
U  consistent  with  having  tlie  2’/lg  below  the  B  ‘is  is  observed  in  the  longer  finite 

polyenes  [3]  is  aliout  .3  eV  2/o).  Our  numerical  results  show  that  thi.'-  ‘crossover" 
point  is  fairly  insensitive  to  system  size  and  to  whctlier  one  works  at  fixed  (h,';  ■  S  or  at 
fixed  l\  .  Note  that  over  much  of  the  range  of  f/  the  energy  of  the  2'.lg  lif's  rougidy  twice 
a.'i  f<i,r  alxive  tiie  ground  state  as  does  the  i  '/A.  slate;  this  strongly  suggests  that  tlie  2‘.4,j 
state  IS  m  some  sense  iomp.is<d  i>!  two  triplets,  as  indicate'd  i>y  the  work  of  i  .ivaii  ami 
.Sehniten  [17], 
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Figure  6  caption;  Electronic  contribution  to  the  energies  of  the  states  I'/lj, 
and  l’i5u  (optical  gap)  as  functions  of  the  Hubbard  U.  Other  parameters  arc  io=2.5eV, 
a  =  4.1eVA,  K  =  51.1eV/A^,  and  Open  chain  boundary  conditions  on  A'=8  sites 

were  used,  and  the  geometry  was  constrained  to  be  uniformly  dimerized.  Note  these 
parameters  give  <5=0.03Afor  (/=4eV,  consistent  with  trans-polyacetylene  For  these  pa 
rameters,  the  falls  below  the  near  f/=5eV. 

B.  Geometries 

To  study  the  geometries  of  the  excited  state';,  we  consider  both  constraint d  structures 
which  are  relevant  to  crystals  of  both  finite  polyenes  and  conjugated  polymers  in  which 
solid  state  effects  impose  external  forces  on  the  individual  chains  [e.g.  octatei  raene,  con¬ 
jugated  polymers  in  solid  stale)  and  re/axcd  structures,  which  are  relevant  to  gas  phase 
experiments  on  finite’  polyenes.  In  Fig.  7  we  show  the  total  energies  .'lee  Ironic  plus 
lattice  elistortion  as  funclicms  of  lh<’  (assumed  constant  !)  dimensionh’ss  eiimerization 
{aSjto)  ferr  tlie  four  slate’s  [dotte’d  in  Fig.  (i.  Two  remarks  concerning  tins  figuie  are 
particularly  significant  hirst,  the  constraint  of  constant  dimerization  (const ramteei  ge 
oiiu’try)  for  the  ex<  iteel  state’s  is  mat  re’ahstic  but  is  nonethelerss  useful  feu  e  onpiansons 
with  fjthe’i  t  lie'oretie  al  re'suits,  siie  fi  as  thoi'-e’ obtanu’d  by  Tavan  and  Schiille'ii  il  7j  Se’ccmd. 
the  a-svinme'lry  be'lvve-en  positive,  anel  ne'gative  value's  erf  the  dimensionie-ss  dime'iization 
arnse’s  simpK  beeause',  toi  llii-.<’  imili' syste>m.'-  with  ojie'ii  boundary  e  erndii  lou'- .  the'  ne'ga 
ti  ve'  vai 111 of  (th !  1 11  e  on  I lo! I . i  to  s v''l  I’ln'.  With  one  fi’wi'r  e ton lile  i lOinis  t  jm 1 1  ; ! a ■  po;-i t.i \‘e 
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values;  these  systems  naturally  have  higher  energies.  1  he  results  shown  in  I'lg.  ^  are  for 
N  =  8  sites,  but  our  results  for  N  =  10-14  show  very  similar  structure,  suggesting  little 
system  size  dependence  in  this  result. 


Figure  7  caption:  Total  energy  as  a  function  of  the  dimensionless  dimerization  aS/io, 
showing  the  shift  in  the  minimum  value  for  various  states  of  importance  in  optical  ex¬ 
periments;  iMj,  B^,2^  Ag,  and  l*i?„  (optical  gap).  Parameters  used  are  U /2t.o=Q.8, 
V/2to=0,  and  /Cto/2Q^=2.0.  Open  chain  boundary  conditions  (b.c.-O)  on  N  —  S  sites 
were  used,  and  the  geometry  was  constrained  to  be  uniformly  dimerized:  see  text 

In  Figs.  8  and  9  we  show  the  relaxed  geometries  on  open  chains  for  the  lowest 
singlet  (Fig.  8)  and  triplet  (Fig.  9)  states.  Note  that  the  abscissa  in  each  of  these  figures 
is  the  number  of  bonds,  so  that  the  total  number  ranges  from  7  (  -  3  <  n  <  3)  to  13 
(~6  S  ”  <  6).  In  Fig.  8,  the  otherwise  uniform  dimerization  shows  a  slight  increase  at  the 
ends;  this  is  consistent  with  data  on  finite  polyenes.  In  contreusl,  Fig.  9  .shows  substantia] 
structure  in  the  dimerization.  In  fact,  in  the  large  system  sizes  the  laMnr  dist.irtion 
is  approaching  that  for  a  soliton  pair.  Thus  it  appears  that  even  at  intei  nn'diale  (e  e) 
coupling,  the  solitonic  nature  of  the  triplet  state  -  or  at  least  of  the  lattice  distortion 
corresponding  to  this  state  remains. 
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Figure  8  caption:  The  relaxed  geometries  of  the  ground  state  of  open  chains  for 
systems  containing  between  seven  and  thirteen  double  bonds,  corresponding  to  finite 
polyenes  from  eight  to  fourteen  {CH)  units  in  the  gas  phase.  Parameters  are  to=2.5eV, 
a  =  4.1eV/A,  K  =  22  -  24eV/A,  and  V=0.  Note  that  K  was  varied  to  eliminate  the 
N  dependence  that  arises  when  the  dimerization  is  constrained  to  be  uniform;  when  one 
allows  relaxation  via  the  small  variation  of  K  values,  the  remaining  system  size  dependence 
is  minimal. 
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Figure  9  caption:  The  relaxed  geometries  of  the  lowest  triplet  state  (13f )  on  open 
chains  for  systems  containing  between  seven  and  thirteen  double  bonds,  corresponding  to 
finite  polyenes  from  eight  to  fourteen  (CH)  units  in  the  gas  phase.  Parameters  are  same 
parameters  as  Fig.  8.  Note  the  lattice  distortion  is  approaching  that  for  a  soliton  pair. 


C.  Optical  Absorption 

In  discussing  the  optical  absorptions,  it  is  particularly  important  to  distinguish  between 
instantaneous  “vertical  transitions”,  either  allowed  or  forbidden,  and  transitions  to  or 
from  states  with  relaxed  geometries,  such  as  those  occurring  in  photo-luminescence.  In 
Fig.  10  we  show  various  energy  gaps  -  labeled  generically  by  “A”  in  the  figure  -  measured 
in  optica!  experiments  as  functions  of  the  Hubbard  U.  BJ'  and  are  used  to 

denote  the  “optically  allowed”  states,  since  for  large  (/,  the  lowest  states  with  non-zero 
optical  absorption  matrix  elements  are  no  longer  the  first  states  of  these  symmetries.  For 
simplicity,  we  have  plotted  on  the  case  with  F  =  0;  hence  there  are  no  excitonic  effects  of 
the  sort  believed  to  be  important  in  poly-diacetylene  cryst.ils  (5). 
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Figure  10  caption  ;  Singlet-singlet  (A5)  and  triplet-triplet  (Ay)  energy  gaps  A 
measured  in  optical  experiments  ais  functions  of  Hubbard  U.  The  parameters  are  <o=2.5eV, 
a  —  4.1eV/A,  K  =  21.0eV/A^,  K=0,  on  an  TV  =  10  site  system,  and  antiperiodic  boundary 
conditions  were  used. 

We  note  that,  as  suggested  by  the  analytic  estimates  in  the  weak  and  strong  coupling 
limits,  the  allowed  triplet-triplet  gap  (Ar)  is  always  less  than  singlet-singlet  gap  (As). 
For  t/  =  0  At  is  almost  exactly  half  As  (as  predicted  by  weak  coupling),  while  for  large 
U,  At  —*  As,  as  expected  in  strong  coupling. 

Finally,  in  Fig.  11  we  plot  the  frequency  dependent  absorption  coefficient  ,  a{u)), 
for  both  singlet-singlet  (bottom  curve,  as(cn))  and  triplet-triplet  (middle  curve,  a'r(‘.<^)) 
absorptions,  as  well  as  their  difference  (top  curve).  The  largest  peak  in  each  case  is  at 
the  edge  of  the  optical  gap,  indicating  that  for  these  intermediate  values  of  U  the  rem¬ 
nant  of  the  one-dimensional  single  particle  density-of-states  peak  survives  in  the  optical 
absorption.  The  prominent  higher  feature  in  the  singlet  absorption  has  been  identified 
previously  [22j  as  a  “decoupled  dimer”  peak:  that  is,  a  feature  that  arises  because  of  the 
dimerization  and  which  increases  in  intensity  with  increasing  dimerization.  This  feature 
is  less  prominent  in  the  triplet  case  because  a  triplet  dimer  has  no  allowed  optical  tran¬ 
sitions,  so  that  the  presence  of  the  single  triplet  on  these  small  systems  tends  to  reduce 
the  intensity  of  this  peak. 
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Figure  11  caption;  The  optical  absorption  coefficients  as  functions  of  frequency  (ui) 
for  allowed  singlet-singlet  transitions  (called  a$  in  the  text),  allowed  triplet-triplet  transi¬ 
tions  (called  ar  in  the  text),  and  for  the  difference  for  an  intermediate  value  of  f/  (=8  eV). 
Other  parameters  are  as  in  Fig.  8.  The  spectra  are  averaged  over  boundary  conditions 
(b.c. =-1,-0. 5,0, 0.5,1)  and  over  system  sizes  A^=8,10,12,14.  The  triplet  and  singlet  absorp¬ 
tions  were  calculated  using  the  relaxed  geometry  of  the  respective  cases.  The  spectra  are 
displaced  for  clarity  and  obtained  from  the  raw  data  by  Lorentzian  broadening;  this  later 
procedure  leaves  some  residual  unphysical  high  frequency  oscillations  in  the  spectra.  Note 
however  the  clear  presence  of  the  decoupled  dimer  peak  in  the  singlet  spectrum. 

V.  Application  to  Photoexcitations  in  Conjugated  Polymers 
A.  The  2Mg  and  states 

In  our  discussion  of  the  energetics  of  excited  states  (cf.  Fig.  6),  we  noted  that  over  much 
of  the  range  of  U  the  energy  of  the  '2}  Ag  lies  roughly  twice  as  far  above  the  ground 
state  as  does  the  1‘^/^u  slate  and  remarked  that  this  strongly  hints  that  the  2Mj  .state 
IS  in  some  sense  composed  of  two  l.nplets,  as  suggested  ahso  by  the  wcmk  of  Tavan  and 
.Schulten  [17].  At  the  sariu’  tirn<\  from  our  studies  of  the  geometries  of  excited  states,  we 
fouml  that  tlie  stalf  had  the  lattice  distortion  appropriate  to  Iwo  neutral  solitons 

fcf.  I‘'ig,  9)  IJoe.s  this  mean  I  hat  one  should  view  tint  slate  a.s  coinjiosed  of  Joxir 
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soiitons?  Such  an  interpretation  was  proposed  in  the  renormalization  grovip  studies  of 
Hayden  and  Mele  [15,  16]  but  is  explicitly  contradicted  in  the  configuration  interactioii 
studies  [17]  and  (at  least  apparently)  in  combined  experimental/theorcticaJ  studies  on 
octatetraene  [24],  Unfortunately,  it  is  easier  to  argue  why  each  of  the  purely  theoretical 
studies  could  be  Incorrect  than  it  is  to  determine  a  definitive  answer  to  this  question. 
In  the  renormalization  group  studies,  the  technique  of  constructing  larger  systems  by 
combining  smaller  ones  produces  a  natural  bias  against  states  with  “long  bonds"  in  the 
valence  bond  description  (see  [18]  for  a  more  thorough  description  of  this  effect).  But 
these  long  bond  states  would  correspond  precisely  to  the  Ag  state  made  up  of  only  two 
triplets,  rather  than  four  soiitons.  On  the  other  hand,  in  the  configuration  interaction 
calculations,  the  geometries  were  not  completely  relaixed,  so  tjhat  one  could  argue  that  a 
bias  existed  against  the  formation  of  soiitons.  At  present,  our  data  on  the  fully  relaxed 
geometry  and  energy  of  the  2Mg  state  on  the  largest  systems  {N  =  12,14)  are  not  sufficient 
to  determine  the  resolution  of  this  question  unambiguously.  This  difficulty  arises  in  part 
because  the  2^Ag  state  is  not  the  lowest  singlet  state  of  Ag  symmetry  and  thus  must  be 
extracted  from  the  Lanezos  data  with  considerable  care. 


B.  Consistency  with  SDPM  data 

Recent  studies  of  the  spin  dependent  photomodulation  (SDPM)  [10,  II]  in  both  degenerate 
ground  state  {trans-{C H)^)  and  non-degenerate  ground  state  polymers  (PPV)  identify 
the  spins  associated  with  long-lived  excitations  and  assign  the  so-called  “high  energy” 
peak  to  triplet-triplet  transitions  associated  respectively  to  neutral  soliton  pairs  in  irans- 
{CH)x  and  to  (spinless)  neutral  bipolarons  (composed  of  two  triplets)  in  PPV.  Although 
more  detailed  quantitative  studies  are  required,  our  present  results  on  triplet  states  and 
their  relations  to  other  excited  states  (such  as  the  2Mj  state)  are  consistent  with  such 
an  interpretation.  If  one  assumes  a  relaxation  mechanism  which  allows  the  optical  ^ 
state  to  cross  over  to  the  lower- lying  2^ Ag  state,  then  one  can  take  advantage  of  the 
schizophrenic  nature  of  this  later  state  to  explain  both  types  of  high  energy  peak.  In  the 
degenerate  ground  state  system  the  “reversed  sense”  ground  state  has  tlie  same  energy, 
so  that  the  long  bond  that  characterizes  the  2^  Ag  state  can  be  weak.  Hence  the  two  spins 
associated  with  this  bond  can  interact  independently  with  the  magnetic  field,  giving  ri.se 
to  a  spin  1/2  signal.  In  non-degenerate  ground  state  polymers,  the  region  spanned  by  the 
long  bond  must  also  span  a  region  of  reversed  sense  of  the  bond  order.  Since  this  i.s  not 
equal  in  energy  to  true  ground  slate,  the  “long- bond"  can’t  be  too  long;  hence  one  .sees  no 
spin  1/2  but  only  spinless  or  triplet  signaturres.  However,  until  substantiated  by  detailed 
quantitative  studies,  sucli  a  theoretical  interpretation  remains  speculatiN  c. 


VI.  Discussion  and  Open  Issues 

We  have  sliown  that  studying  I’eierls  <>\teiided  Hubbard  models  with  Lain /os  exai  t  diag 
onalizatioii  t('rhnH)ue.s  and  bouridarr'  <  omlilion  .tveiaging  metlioris,  allows  one  to  explore 
the  pifiperties  of  triplet  stales  m  both  fimtf'  jiolvene.s  and  conjugated  poKuioiS  over  tli<' 
entile  range  of  p  and  e  <■  < ,  Miphng.s,  Om  re.'-iill,--  are  <  onsistent  wit.li  the  kn'>wn  .nialytic 
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limits.  More  extensive  studies  should  provide  detailed  interpretations  of  data  from  real 
finite  polyenes,  both  in  gas  and  crystal  phases,  and  also  from  conjugated  polymers. 

Apart  from  these  more  detailed  studies,  there  are  several  additional  open  issues  which 
we  plan  to  address.  First,  we  shall  attempt  to  find  a  definitive  resolution  of  the  problem 
of  the  geometry  of  the  2Mp  state.  Second,  we  intend  to  investigate  quantitatively  the 
consequences  of  “confinement”  -ie,  the  lack  of  ground  state  degeneracy  -  in  conjugated 
polymers.  Third,  we  propose  to  study  the  potential  roles  of  possible  charged  triplet  states 
in  doped  conjugated  polymers.  Finally,  we  plan  to  explore  the  manifold  of  higher  spin 
magnetic  states  in  these  models  and  systems,  in  the  hope  that  they  may  provide  still 
further  insight. 
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